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ABSTRACT 


f\  ^  |  .  .  .1.  I 

7-fn  this  paper  we  determin&the  nonnegativity  of  the  principal  components  of  an  n  x  n 
nonnegative  matrix  P  in  terms  of  the  marked  reduced  graph  %{A)  of  A  =  P  —  the 

minus  Mfmatrix  which  can  be  associated  v/ith  P.  We  then  apply  this  result  to  consider  var¬ 
ious  types  of  nonnegative  bases  for  the  Perron  eigenspace  of  P  which  can  be  extracted  from 
a  certain  nonnegative  matrix  which  is  a  polynomial  in  P.  We  also  obtain  a  characterization 
for  the  eigenprojection  on  the  Perron  eigenspace  of  P  to  be,  itself,  a  nonnegative  matrix. 
Our  results  provide  new  proofs  and  extensions  of  results  of  Friedland  and  Schneider  and 
of  Hartwig,  Neumann,  and  Rose. 


AMS  (MOS)  Subject  Classifications:  15A18,  15A21,  14A48,  05C28 

Key  Words:  minus  M-matrix,  nonnegative  matrix,  principal  component,  transform, 
component,  eigenprojection,  Perron  eigenspace,  nonnegative  basis, 
semi-preferred  basis,  proper  combinatorial  basis,  asymptotic  behavior 
.of  powers 
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1  INTRODUCTION 


In  this  paper  we  determine  the  nonnegativity  structure  of  the  principle  com¬ 
ponents  of  an  n  x  n  nonnegative  matrix  P  in  terms  of  the  marked  reduced 
graph  72(i4)  of  the  associated  minus  M-matrix  A  =  P  -  p(P)I<  We  then  ap¬ 
ply  this  result  to  obtain  a  characterization  for  when  the  eigenprojection  on 
the  generalized  eigenspace  corresponding  to  the  Perron  root  (subsequently 
called  the  Perron  eigenspace  of  P  or  A)  is  itself  a  nonnegative  matrix. 
We  also  apply  this  result  to  obtain  new  proofs  and  extensions  of  results 
of  Hartwig,  Neumann,  and  Rose  [3]  and  Friedland  and  Schneider  [5] 

It  is  well  known  that  the  Perron  eigenspace  of  P  has  a  basis  of  nonneg¬ 
ative  vectors  with  many  specified  properties,  see  Rothblum  [11],  Richman 
and  Schneider  [10],  Schneider  [12],  and  Hershkowitz  and  Schneider  [6],  The 
approach  taken  in  these  papers  is  matrix  combinatorial  and  uses  the  Frobe- 
nius  normal  form.  More  recently  the  authors  of  [3]  gave  an  analytic  proof 
for  the  existence  of  a  nonnegative  basis.  They  found  a  nonnegative  matrix 
J  which  turns  out  to  be  a  polynomial  in  P  and  whose  columns  contain  a 
basis  for  the  Perron  eigenspace  of  P.  This  suggests  that  the  two  approaches 
are  related.  In  this  paper  the  relation  is  investigated. 


We  now  describe  the  contents  of  the  paper  in  more  detail.  Section  2  con¬ 
tains  most  of  the  notions  which  we  use  in  the  paper.  In  Section  3  we  prove 
our  first  main  result  (Theorem  1):  Let  P  be  an  nxn  nonnegative  matrix  in 
Frobenius  normal  form  and  let  £(*)  be  the  k-th  principal  component.  Let  i 
and  j  be  vertices  in  TZ(A)  and  let  d  =  d(i,j )  be  the  (singular)  distance  from 
i  to  j.  Then  the  ( i,j)-th  block  of  £(*)  is  strictly  positive  if  d  =  k  +  1.  We 
mention  that  according  to  Lemma  2(i),  the  (i,  j)-th  block  of  #(*)  is  zero  if 
d  <  k  +  1.  The  results  of  [3]  are  a  corollary  of  this  result.  Another  outcome 
of  this  result  is  that  a  nonnegative  basis  can  be  extracted  from  the  columns 
of  J  which  is  strongly  combinatorial  in  the  sense  defined  in  Section  2. 


In  Section  4  we  characterize  when  the  eigenprojection  Z  =  Z(°)  on  the 
Perron  eigenspace  of  P  is  itself  a  nonnegative  matrix.  We  begin  by  showing 
that  a  necessary  condition  for  this  to  happen  is  that  in  TZ(A)  no  nonsingular 
vertex  lies  on  the  interior  of  a  simple  path  connecting  two  singular  vertices. 
This  leads  to  our  characterization  (cf.  Theorem  2)  for  Z  to  be  nonnega¬ 
tive  in  terms  of  the  existence  of  special  nonnegative  bases  for  the  Perron 
eigenspaces  of  A  and  AT  which  are  bi-orthogonal  and  whose  distinguishing 
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feature,  compared  with  other  types  of  nonnegative  bases  which  the  Perron 
eigenspaces  of  these  matrices  possess,  is  their  comparative  sparseness. 

In  Section  5,  using  Theorem  1,  we  re-prove  results  contained  in  [5]  on 
the  asymptotic  behavior  of  the  powers  of  P.  In  the  special  case  when  all 
the  vertices  in  P(A)  come  from  primitive  diagonal  blocks  in  the  Frobenius 
normal  form  of  F,  'he  result  asserts  that  (Pm);j  behaves  asymptotically  as 
7 mdsmziS~1\P{i. .?}),  where  7  is  a  positive  constant,  6  =  S(i,j)  is  the  local 
distance  between  from  i  to  j,  s  is  the  local  spectral  radius  of  P{i,j},  the 
principal  submatrix  of  P  index  by  all  vertices  /  which  lie  on  a  path  from 
i  to  j  in  7 £(P),  and  7iiS~ ^{P{i,j})  is  the  (£  -  l)-th  principal  component 
of  P{i,j }  corresponding  to  the  eigenvalue  s.  In  the  general  case,  without 
any  primitivity  assumption,  wo  consider  the  sequence  MPm ,  where  M  is  a 
smoothing  matrix  which  is  a  polynomial  in  P. 

In  this  introduction  we  have  described  our  results  in  terms  of  the  non¬ 
negative  matrix  P.  In  Sections  2-4  it  will  be  convenient  to  state  our  result 
in  terms  of  the  associated  minus  M-matrix  A  given  above. 


2  NOTATIONS  AND  PRELIMINARIES 

For  a  positive  integer  n  we  denote  by  (n)  the  set  {1, . . n}. 

For  an  nx  n  matrix  A  we  denote  by: 

N(A)  -  the  nullspace  of  A. 

E(A)  -  the  generalized  nullspace  of  A,  viz.  ^ {(An ). 

u(A)  -  the  index  of  0  as  an  eigenvalue  of  A,  viz.  the  size  of  the  largest 
Jordan  block  associated  with  0. 

Let  a  C  (n).  By  A[a]  we  shall  denote  the  principal  submatrix  of  A 
whose  rows  and  columns  are  determined  by  or.  Similarly,  for  an  n-vector 
x,  we  shall  denote  by  z[a]  the  subvector  of  x  whose  entries  are  indexed  by 
a.  For  an  array  C  we  shall  use  C  >  0  to  denote  when  all  its  entries  are 
nonnegative  numbers.  C  >  0  shall  denote  the  fact  that  C  >  0,  but  C  ^  0. 
C  >  0  shall  denote  the  fact  that  all  of  the  entries  of  C  are  positive  numbers. 


In  all  our  considerations  we  shall  assume  that  A  is  an  n  x  n  real  matrix 
given  in  a  block  lower  triangular  form  with  p  square  diagonal  blocks  as 
follows 


(A i,i  0  ...  0  \ 

I  ^2,1  -^2,2  0  j 


(2.1) 


\  Ap%  i  .  ' 


where  each  diagonal  block  is  irreducible  or  the  lxl  null  matrix.  The  above 
form  is  the  so  called  Frobenius  normal  form  of  A.  It  is  well  known  that 
any  square  matrix  is  symmetrically  permutable  to  such  a  form.  The  reduced 
graph  of  A,  TZ(A),  is  defined  to  be  the  graph  with  vertices  l,...,p,  where 
(t,j)  is  an  arc  from  i  to  j  if  Aij  £  0.  A  vertex  i  in  TZ(A)  is  said  to  be 
singular  if  A{ti  is  singular.  Otherwise  the  vertex  is  called  nonsingular.  The 
set  of  all  singular  vertices  in  7£(/l)  will  be  denoted  by  <5(/l).  A  sequence  of 
vertices  (ij, . . . ,  4)  in  FL(A)  is  said  to  be  a  path  from  4  to  4  if  there  is  an 
arc  in  'R.(A)  from  ij  to  ij+i,  Vj  6  {k  -  1).  The  path  is  said  to  be  simple 
if  4, . .  .,4  are  distinct.  The  empty  path  will  be  considered  a  simple  path 
linking  every  vertex  ?  €  F-(A)  to  itself.  If  there  is  a  path  (in  7 £(A))  from  i 
to  j  we  shall  write  that  i  y  j.  If  i  £  j  and  there  is  a  path  from  i  to  j  we 
shall  write  that  i  y  j.  If  there  is  a  path  from  i  to  j  define  the  ( singular ) 
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distance  d(i,j)  from  i  to  j  to  be  the  maximal  number  of  singular  vertices 
on  a  simple  path  connecting  i  and  j.  If  there  is  no  path  from  i  to  j  we  set 
d(i,j )  =  —1.  In  particular  it  follows  from  our  definition  that  d(i,  i)  =  0  if  i 
is  a  nonsingular  vertex  and  d(i,  t)  =  1  if  i  is  a  singular  one.  Note  that  each 
vertex  i  6  7£(A)  can  be  thought  of  a  subset  of  (n)  which  consists  of  those 
elements  in  (n)  upon  which  the  i-th  diagonal  block  in  A  is  indexed. 

If  0  is  an  eigenvalue  of  A,  then  in  a  punctured  neighborhood  of  0  which 
contains  no  other  eigenvalues  of  A,  the  resolvent  operator  (e/  -  A)~ 1  admits 
the  Laurent  expansion 


(d-A)-'  =  £  =$■  +T(€),  (2.2) 

k= 0 

where  Zt\  is  the  eigenprojection  of  A  on  E(A)  (that  is,  the  projection 
on  E(A)  along  the  join  of  all  eigenspaces  of  A  corresponding  to  eigenval¬ 
ues  other  than  zero)  and  where  T(e)  is  an  analytic  operator  in  e  defined 
throughout  the  nonpunctured  neighborhood  of  0  satisfying  T(0)  =  0  and 
ZaT(c)  =  T(c)Zy\  =  0  (cf.  Kato  [8,  pp. 34-43]).  We  call  the  matrices 
Z(k)  :=  A^'Za,  k  =  0, 1, . ..,!/( A)  —  1,  the  principal  components  of  A  ( cor¬ 
responding  to  the  eigenvalue  0).  Recall,  e.g.  Lancaster  and  Tismenetsky 
[9],  that  all  principal  components  of  A  are  functions  of  A  and  hence  also 
polynomials  in  A.  Note,  however,  that  the  principal  components  as  defined 
here  differ  by  factorial  multiples  from  those  introduced  in  Lancaster  and 
Tismenetsky  (see  [9,  p.314]). 

Consider  Tl(A).  We  shall  let 

(i,j)  =  {/  €  TZ(A)  :  it  It  j} 

and  we  shall  let  A {i,j}  =  /l[(i,y)].  The  following  lemma  is  readily  obtained: 

LEMMA  1  For  every  polynomial  D  =  q(A)  we  have  that  D[(i,j)]  =  g(A{i,j}). 
In  particular,  as  all  the  principal  components  of  A  are  polynomials  in  A, 
ZW((i,i»  =  ((/!){.•, 3))%,U). 


Proof:  Partition  B  =  q(A)  in  conformity  with  A  and  compute  B[(i,j)] 
from  the  powers  of  A  which  appear  in  q(A).  □ 
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Let  P  be  an  n  x  n  nonnegative  matrix.  The  Perron  Frobenius  theory 
(cf.  Berman  and  Plemmons  [l])  tells  us  that  the  spectral  radius  of  P,  given 
by  the  quantity 


p(P)  =  max{|A|  :  det(P  -  AJ)  =  0} 

is  an  eigenvalue  of  P  to  which  there  corresponds  a  nonnegative  eigenvector. 
In  particular,  if  P  is  irreducible,  then  p(P)  is  simple  and  the  corresponding 
eigenvector  is,  up  to  a  multiple  by  a  scalar,  positive.  The  matrix  A  — 
P  —  p(P)J,  which  has  all  its  off-diagonal  entries  nonnegative,  is  the  n  x  to 
minus  M-matrix  which  we  associate  with  P  and,  in  several  sections  of  our 
paper,  it  will  be  convenient  to  work  with  A  rather  than  with  P.  (We  call  A  a 
minus  M-matrix  if  -A  is  an  M-matrix.  For  the  many  equivalent  conditions 
for  a  real  matrix  with  nonpositive  off-diagonal  entries  to  be  an  M-matrix 
see  Berman  and  Plemmons  [l,  Chp.6].)  Suppose  now  that  m  =  dim(E(A)). 
It  is  known  that  m  is  equal  to  the  number  of  singular  vertices  in  7Z(A),  e.g., 
Cooper  [2].  Rothblum  [11]  has  shown  that  u(A)  is  equal  to  the  maximum 
over  all  lengths  of  the  simple  paths  in  TZ(A),  a  result  to  which  we  shall 
refer  as  the  Roiliblum  index  theorem.  Let  <S(-4)  =  {c*i , . . .  ,am}.  Rothblum 
[11]  and,  independently,  Richman  and  Schneider  [10]  have  shown  that  E(A) 
possess  a  basis  of  nonnegative  vectors  having  the  following 

properties: 

uW[i\  = 

We  shall  call  a  nonnegative  basis  for  E(A)  satisfying  (2.3)  strongly  combi¬ 
natorial. 


»  0, 
0, 


iff  i  X  aj, 
otherwise. 


(2.3) 


Let  A  be  a  minus  M-matrix  given  in  the  form  (2.1).  Then  AT  is  now 
a  minus  M-matrix  in  upper  triangular  Frobenius  normal  form.  If  we  were 
to  introduce  formally  the  concept  of  the  reduced  graph  for  matrices  in  up¬ 
per  triangular  form,  we  would  observe  that  in  going  from  A  to  AT,  the 
direction  of  the  access  between  the  vertices  is  reversed.  In  any  event,  by 
possibly  transforming  AT  to  a  lower  triangular  form  via  similarity  permu¬ 
tations,  it  follows,  on  applying  the  aforementioned  results  for  nonnegative 
bases  for  matrices  in  lower  triangular  Frobenius  normal  form,  that  E(AT) 
has  a  nonnegative  basis  of  vectors  such  that  for  i  e  ( p )  and 


j  G  («»), 


iff  in  7?.(A),  aj  >  i, 
otherwise. 
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Motivated  by  this  observation  and  by  Hershkowitz  and  Schneider  [6]  we  now 
introduce  the  notions  of  column  and  row  proper  combinatorial  bases. 

DEFINITION  1  Let  A  be  the  n  X  n  minus  M-matrix  given  in  (2.1)  and 
consider  11(A).  A  nonnegative  basis  . . . ,  is  called  a  proper  column 
combinatorial  basis  for  E(A)  if 

>0  =*>  ibotj 


and 

>  0 

for  all  i  G  (p)  and  j  G  (m).  Similarly,  a  nonnegative  basis  vW, ..  for 

E(At)  is  rat  ltd  v  proper  row  combinatorial  basis  for  E(AT)  if 

>  0  =>  aj  yi 

and 

u(i)[aj]  »  0, 

for  all  i  G  (p)  and  j  G  (m). 

We  end  this  section  with  the  following  observation  which  will  be  useful 
for  our  results  in  the  next  sections: 

LEMMA  2  Let  A  be  as  in  (2.1)  and  let  i  and  j  be  vertices  in  11(A). 

(i)  For  all  k  >  d(i,j),  Z$  =  0. 

(ii)  If,  in  addition,  A  is  a  minus  singular  M-matrix,  then  > 

0. 


Proof:  (i)  By  the  result  in  Friedland  and  Hershkowitz  [4]  and  Her¬ 
shkowitz,  Rothblum,  and  Schneider  [7]  the  index  of  A{i,j)  does  not  exceed 
k.  Hence,  by  Lemma  1,  Z^[{i,j)]  =  0  and  the  result  follows. 

(ii)  This  is  an  consequence  from  the  resolvent  expansion  of  i4{i,j}  in 
a  sufficiently  small  punctured  neighborhood  of  0  considering  that  the  fact 
that  (el  -  i4{i,y})-1  >  0,  Ve  >  0.  □ 
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3  PRINCIPAL  COMPONENTS 


We  mentioned  in  the  introduction  that  the  first  approach  to  the  questions  of 
existence  of  a  nonnegative  basis  to  the  Perron  eigenspace  of  a  nonnegative 
matrix  P  and  the  properties  of  such  a  basis  were  matrix  combinatorial  in  the 
sense  that  they  were  based  on  the  Frobenius  normal  form.  This  approach 
was  developed  by  Rothblum  [11]  and  Richman  and  Schneider  [10],  A  later 
approach  to  these  questions,  obtained  by  Hartwig,  Neumann,  and  Rose  [3], 
is  analytic  in  the  sense  that  it  utilizes  the  resolvent  expansion,  but  does  not 
involve  the  Frobenius  normal  form.  In  [3]  the  authors  show  that  if  Zp_p(p)j 
is  the  eigenprojcction  on  the  Perron  eigenspace  of  P,  then  for  A  —  p(P)  >  0 
sufficiently  small,  the  matrix 

J  =  (AJ-Pr^pj;  (3.1) 

is  nonnegative  and  its  columns  contain  a  basis  for  the  Perron  eigenspace.  We 
now  observe  that  if  P  is  in  lower  block  triangular  Frobenius  normal  form, 
then  since  J  is  a  polynomial  in  P,  J  is  also  a  block  lower  triangular  matrix 
whose  partitioning  conforms  with  that  of  P.  This  suggests  that  there  is  a 
connection  between  the  combinatorial  and  the  analytic  approaches  to  the 
existence  of  a  nonnegative  basis  described  above.  In  this  section  we  shall 
make  the  connection  more  precise.  It  will  be  convenient  for  us  to  work  and 
state  our  results  in  terms  of  the  associated  minus  M-matrix  A  =  P  —  p{P)I. 

Our  first  main  result  is  the  following: 

THEOREM  1  Let  A  be  a  minus  M-matrix  in  Frobenius  normal  form  given 
by  (2.1)  and  set  Z  =  Za ■  Suppose  i,j  £  P(A).  If  d(i,j)  =  k  >  1,  then 

zg_,)  >  0. 

Proof:  The  theorem  is  proved  by  induction  on  number  of  blocks  p  in 
the  Frobenius  normal  form. 

Let  p  =  1.  Then  i  =  j  =  1.  If  1  is  a  nonsingular  node,  then  rf(l,  1)  =  0 
and  hence  of  no  interest.  If  1  is  a  singular  node,  then  d(l ,  1)  =  1  and,  as  is 
well  known,  Z^0)  >  0.  This  proves  the  case  p  =  1. 

Assume  now  that  the  result  is  true  for  all  minus  M-matrices  with  fewer 
than  p  blocks  in  their  Frobenius  normal  form  and  suppose  that  A  is  a  minus 
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M-rnatrix  with  p  blocks  in  its  normal  form.  Let  ttj  6  7S(/1)  be  such  that 
d(i,j)  —  k  >1.  We  need  to  distinguish  four  cases. 

Case  1.  i  and  j  are  nonsingular  vertices. 


(a)  Since  d(i,j)  =  k  there  exists  a  simple  path  from  i  to  j  containing  k 
singular  vertices.  Let  (i,  m)  be  the  first  arc  in  this  path  so  that  Ai>m  >  0. 
Then  d(m,j)  =  k.  Thus,  by  our  inductive  assumption  and  Lemma  2  applied 
to  A{m,  j},  we  have  that  Z^f1^  >  0.  Furthermore,  again  by  Lemma2,  since 
the  index  of  A{i,j)  is  k ,  Z^k~^[{i,j)]  >  0  and  ZW[{i,j)]  =  0.  Now  as 


and  as  every  term  in  this  sum  except  for  the  last  one  is  nonnegative  with  at 
least  one  of  the  terms  being  semipositive,  it  follows  that 


Rut  —A{t{  is  a  nonsingular  and  irreducible  M-matrix  and  so  its  inverse  is 

strictly  positive.  It  follows  that  Zjj~^  >  0  and  that  each  of  its  columns  is 
either  strictly  positive  or  zero. 


(b)  Let  (l,j)  be  the  last  arc  in  a  path  from  i  to  j  which  contains  k 
singular  vertices.  As  in  part  (a)  we  obtain  that  Aij  >  0  and  it  follows  by 
our  inductive  assumption  that  zfj~^  >  0.  As  before,  from 


0  =  zg>  =  Zfr^Aij  + . . .  +  zg-"/ 1«  +  . . .  +  zg-’^u,  (3.4) 

it  follows  that 

(3.5) 


By  an  argument  similar  to  the  one  following  (1.3)  it  follows  U'.t  zjj~ ^  >  0 
and  that  each  of  its  rows  is  either  zero  or  positive. 


Combining  the  results  of  (a)  and  (b)  we  obtain  that  >  0. 

Case  2.  i  :  a  nonsingular  vertex  and  j  is  a  singular  one. 


(a)  As  in  Case  1  we  choose  a  simple  path  from  i  to  j  containing  k  singular 
vertices.  By  the  argument  of  Case  1(a),  which  did  not  use  the  nonsingular¬ 
ity  of  the  vertex  j  we  prove  that  z\^  >  0  and  each  of  its  columns  is  either 
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positive  or  zero. 


(b)  Choose  /  as  in  Case  1(b).  Since  d(l ,  i)  =  A:  —  1  it  follows  from 
our  inductive  assumptions  and  Lemma  2  that  Z(fc-1)[(/,i)]  =  0  and  that 
Z-j  =  0.  Hence 

0  =  z\j  =  (3-6) 

Thus  each  row  of  Z is  a  left  null  vector  of  the  irreducible  M-matrix 
-Ajj  and  hence  is  either  positive  or  zero. 

Combining  the  results  of  (a)  and  (b)  we  obtain  that  >  0. 

Case  3.  i  is  a  singular  vertex  and  j  is  a  nonsingular  one. 


The  proof  of  this  case  is  similar  to  the  proof  of  Case  2  with  the  roles  of 
the  vertices  i  and  j  reversed. 

Case  4.  i  and  j  are  singular  vertices. 


Observe  that  by  Lemma  2  the  trailing  principal  submatrix  of  2f(*-1)[(i,  j)] 
obtained  by  deleting  the  first  block  row  and  the  first  block  column  of 
Z(fc-1)[(i,  j)]  is  zero.  A  similar  argument  shows  that  the  leading  principal 
submatrix  of  j)]  obtained  by  deleting  the  last  block  row  and  the 

last  block  column  of  Z(fc-1)[(i,  j)]  is  zero.  Hence  all  the  blocks  of  #(*-1)[(i,j)] 
with  the  possible  exception  of  Z-j~^  are  zero.  However,  Z^k~^[{i,j)}  ^  0, 
so  that  by  Lemma  2(ii)  we  must  have  that  zjj~^  >  0.  Repeating  the  proof 

of  Case  2(b)  we  obtain  that  each  row  of  Z-j~^  is  either  positive  or  zero. 
Similarly,  by  reversing  .he  roles  of  i  and  j  in  this  proof,  we  obtain  that  each 
column  of  zjj is  either  positive  or  zero.  Together  these  conclusions  give 
that  >  0.  □. 


Theorem  1  has  the  immediate  consequence  that  on  choosing  e  =  A  - 
p(P )  >  0  sufficiently  small,  the  nonnegative  basis  for  E(A),  A  =  P  -  p(P)I, 
which  can  be  extracted  from  the  columns  of  the  matrix  J  in  (3.1)  is  a 
combinatorial  basis.  To  see  this  note  that  from  (2.2)  and  (3.1)  it  follows 
that 


J  - 


i<A)- 1 

E 


k—0 


Z(V 


(3.7) 
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and  recall  that  as  A  is  in  lower  triangular  Frobenius  normal  form,  J  itself 
is  lower  triangular.  The  nonzero  diagonal  blocks  of  J  are  positive  and  oc¬ 
cur,  precisely,  ir.  position^  corresponding  to  the  singular  blocks  of  A.  The 
remaining  diagonal  blocks  of  J  are  zero.  Now  from  each  block  column  of 
J  containing  a  positive  diagonal  block  of  J  select  a  column  of  J.  Clearly 
the  m  columns  thus  selected  are  linearly  independent  and  form  a  basis  for 
E(A)  which,  by  Theorem  1  and  simple  asymptotic  consideration,  satisfies 
the  requirements  of  (2.3)  for  being  a  combinatorial  basis. 


Continuing,  we  mention  that  both  Rothblum  [11]  and  Richman  and 
Schneider  [10]  show  that  nonnegative  bases  for  E(A)  can  be  selected  possess¬ 
ing  stronger  properties  than  just  (2.3).  In  particular,  Richman  and  Schneider 
show  that  E(A)  has  a  so  called  preferred  basis  which  is  a  set  of  nonnegative 
vectors  ul1), . . . ,  which  forms  a  basis  for  E(A)  satisfying  (2.3)  with  the 
added  stipulation  that 

m 

AuW  = 

i=o 

where  the  c*(j’s  are  scalars  satisfying 


f  >  0,  if  aj  s-  a*, 
\  0,  otherwise. 


(3.8) 


Thus,  in  view  of  the  foregoing  paragraph,  it  should  be  noted  that  one  can 
find  examples  such  that,  even  for  e  >  0  sufficiently  small  in  (3.7),  a  basis 
extracted  from  the  columns  of  J  cannot  be  a  preferred  basis  as  the  following 
example  shows:  Let 

/0  0  ON 

A  =  (  1  0  0  . 

\0  1  0/ 

As  dim(£(A))  =  3,  we  know  that  the  three  columns  of  J  form  a  basis  for 
E(A).  Let  J t'l,  i  =  1,2,3,  denote  the  columns  of  J.  Now  since  AJ  -■  JA, 
we  see  that 

•AjW  =  c,,tJ«  +  Cl, 2*^  +  cl,3>^^ 


with 


Ci,i  =  0,  Ci,2  =  1,  and  cli3  -  0. 


We  observe  therefore  that,  while  3  y  1,  ci,3  ^  0  defying  the  requirements 
of  (3.8).  However,  we  can  show  that,  for  sufficiently  small  e  >  0,  a  basis 
satisfying  a  somewhat  weaker  conditions  than  being  preferred  can  be  chosen. 
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DEFINITION  2  Let  A  be  an  n  y  n  minus  M-matrix  given  in  the  form 
(2.1).  A  nonnegative  basis  . . . ,  for  E(A)  is  said  to  be  semi-preferred 

if  it  satisfies  the  requirements  of  (2.3)  and  if  the  following  conditions  hold: 

m 

Au W  =  =  (3.9) 

i= i 

where  cj.j  >  0,  k,j  =  Further  we  have  the  implication: 


k  ^  j  and  Ck,j  >  0  aj  >~  a;..  (3.10) 


We  now  have  the  following  corollary  to  Theorem  1: 

COROLLARY  1  Suppose  that  A  is  a  minus  M-matrix  given  in  the  form 
(2.1).  Then  for  sufficiently  small  e  >  0  a  basis  can  be  extracted  from  the 
columns  of  J  given  in  (3.7)  which  is  semi-preferred. 

Proof:  From  (2.2),  the  definition  of  the  Z^'s  in  Section  2,  and  (3.7)  we 
see  that  e  >  0  sufficiently  small  can  be  chosen  so  that 


"('t)-l  y(k) 

"  =  £  V  - 

j=l 


As  before  let  be  the  singular  vertices  in  7Z(A).  For  each  j  = 

1  choose  from  the  aj-th  block  column  of  J  a  column  of  J,  say  it 

is  column  fij  of  J.  Denote  the  columns  of  J  so  chosen  by 
Suppose  now  that  Ck,i,...,Ck,m  are  nonnegative  numbers  such  that 


AJM  = 

j= i 

Let  q  ^  k  and  suppose  that  ck,q  >  0.  Then  we  can  write  that 
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0  < 

=  Z%,  «w>'W 


_  v-'m 

-  2-,j=l 


<*jl io1-1 


-  2^s= o 


+ Eio1-1 


+E;:!«JE"io)‘1 


Let  j  G  (<7  —  1)  and  consider  the  inner  summation  sign  in  the  second  term 
of  the  expression  immediately  above.  If  otj  is  a  vertex  in  TZ(A)  such  that 
aq  )f-  ajy  then  all  summands  under  this  summation  are  0  because  Za^Qj  =  0 
for  all  5  ss  0, ... ,  v(A)  -  1.  On  the  other  hand,  if  aq  >~  aj,  then  from  The¬ 
orem  1  we  can  deduce  that  the  nonzero  coefficient  of  the  lowest  exponent 
in  e  >  0  is  a  positive  vector  for  it  is  a  column  of  a  positive  submatrix  of  an 
appropriate  principal  component  of  A.  Hence,  for  each  such  j  we  can  choose 

c  >  0  sufficiently  small  so  that  >  0-  Overall,  we  can 

choose  e  >  0  sufficiently  small  so  that  the  entire  second  term  in  the  above 
expression  is  nonnegative.  As  Z^q\aq  >  0,  (2f(°))Kl[a9]  >  0,  showing,  in 
turn,  that  for  some  s  >  1,  >  0.  Hence  d(aq,ak)  ^  0  so  that  aq  has 

access  to  a*.  This  completes  the  proof.  □ 


We  finally  comment  that  Theorem  1  and  Corollary  1  strengthen  the 
results  of  Theorem  1  of  Hartwig,  Neumann,  and  Rose  [3]  which  were  briefly 
described  in  the  beginning  of  the  section.  There  the  authors  showed  that 
a  nonnegative  basis  for  E(A)  could  be  extracted  from  the  columns  of  J. 
Theorem  1  and  Corollary  1  above  now  show  that  such  a  basis  can  be  chosen 
which  possesses  a  variety  of  combinatorial  properties. 
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4  NONNEGATIVE  EIGENPROJECTIONS 


In  this  section  we  shall  develop  necessary  and  sufficient  conditions  for  the 
eigenprojection  on  the  Perron  eigenspace  of  an  n  x  n  nonnegative  matrix  P 
to  be,  itself,  a  nonnegative  matrix.  We  shall  see  that  the  nonnegativity  of 
the  eigenprojection  is  tantamount  to  the  existence  of  certain  “sparse”  non¬ 
negative  bases  for  the  Perron  eigenspaces  of  P  and  PT,  respectively.  We 
mention  that  the  appearance  of  the  Perron  eigenspace  of  PT  as  well  as  that 
of  P’s  is  of  no  surprise  since  Za  >  0  if  and  only  if  Z^t  =  (Za)t  >  0.  Again 
it  will  be  convenient  to  state  and  prove  our  results  in  terms  of  the  associated 
minus  M-matrix  A  —  P  —  p(P)I. 

To  begin  with  recall  that  A  is  in  lower  triangular  Frobenius  normal 
form  with  p  diagonal  blocks.  Consider  71(A).  By  Lemma  1,  for  any  vertex 
i  e  7 2(A),  Za[i\  =  Za{{}  —  Za[{]-  Thus,  as  generally  known,  but  confirmed 
by  our  Theorem  1, 

7  pi  —  /  ^  0,  if  i  is  a  singular  vertex  in  72(A), 

A  *  1 0,  otherwise. 

It  follows  that  the  matrix  Za  -  I  is  a  minus  M-matrix  in  (lower  triangular) 
Frobenius  normal  form  with  q  >  p  diagonal  blocks.  The  singular  diagonal 
blocks  oiZA-I  are  of  same  size  and  occur  in  same  positions  as  the  singular 
diagonal  blocks  of  A.  The  nonsingular  diagonal  blocks  of  Za  —  T,  which  are 
now  all  1  x  1,  occur  at  positions  which  correspond  to  diagonal  entries  within 
nonsingular  diagonal  blocks  of  A.  The  implications  of  these  facts  can  be 
expressed  relative  to  the  reduced  graphs  of  A  and  Za  —  /•  For  that  purpose, 
as  before,  denote  by  ai , . . . ,  am  the  singular  vertices  of  7l( A)  and  now  in¬ 
troduce  Piy..-, Pm  to  denote  the  singular  vertices  of  1Z(Za  -  /)•  Then  as 
subsets  of  indices  from  (n),  a,-  =  i  —  1  ,...,m.  Furthermore,  as  subsets 
of  indices  from  (n),  every  nonsingular  vertex  of  71(Za  —  I)  is  contained  in 
some  nonsingular  vertex  of  7£(A). 

Next  we  find  a  simple  lemma  which  gives  a  necessary  condition  for  Za 
to  be  nonnegative. 

LEMMA  3  Let  A  be  an  nxn  minus  M-malrix.  Then  a  necessary  condition 
for  Z  =  Za  to  be  nonnegative  is  that  in  71(A)  no  nonsingular  vertex  lies  on 
the  interior  of  a  simple  path  connecting  two  singular  vertices. 
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Proof:  Suppose  that  the  statement  is  false.  Then  there  exist  vertices 
i,j  £  5(A)  and  k  g  <S(A)  such  that 

i  y  k  >-  j 


and  such  that 

d(i,k)  =  d(k,j)  =  1. 

But  then,  by  Theorem  1, 

Zitk  >  0  and  Zk,j  >0. 

This  means,  according  to  the  expla  .ation  preceding  the  lemma,  that  in 
7 Z(Z  - 1)  there  is  a  path  of  length  at  least  2,  contradicting  the  fact  that  Z 
is  a  nonnegative  matrix  whose  Perron  root  has  index  1.  □ 

We  are  now  ready  to  state  the  main  result  of  this  section. 


THEOREM  2  Let  A  be  an  n  x  n  minus  M-matrix.  Then  necessary  and 
sufficient  conditions  for  Z  =  Zt\  to  be  nonnegative  are: 

(i)  In  7 £(A)  no  nonsingular  vertex  lies  on  the  interior  of  a  simple  path 

connecting  two  singular  vertices. 

(ii)  E(A)  and  E(AT )  have  nonnegative  column  and  row  combinatorial  bases 

and  respectively,  which,  fori  £  (p)  and 

j  6  (m),  possess  the  following  properties: 

(a)  u(i)[i]  >  0  =»  i  =  otj  or  i  is  a  nonsingular  vertex. 

(b)  >  0  =i>  i  =  ctj  or  i  is  a  nonsingular  vertex. 

Proof:  We  begin  by  proving  the  necessity  part.  Suppose  therefore  that 
Z  >  0.  Then  (i)  holds  by  Lemma  3.  Consider  now  the  minus  M-matrix 
Z-I.  Because  Z  is  a  projection,  we  know  that  v(Z-I)  <  1.  The  Rothblum 
index  theorem  now  implies  that  the  length  of  any  simple  path  in  7 Z(Z-I)  is 
at  most  1.  Hence  for  any  singular  vertex  aj  6  TZ(Z-I),  a  vertex  i  £  Tl(Z—I) 
satisfies  that  d(i,otj )  >  0  if  and  only  if  i  =  otj  or  i  is  a  nonsingular  vertex. 
Thus,  by  (2.3),  E(Z-I)  has  a  nonnegative  proper  column  combinatorial  ba¬ 
sis  u(a), . . . ,  satisfying  the  requirements  in  (a)  with  respect  to  7 Z(Z  —  I). 
From  E(A)  =  E(Z  -  I)  and  from  the  relation  between  7 Z(Z  - 1)  and  11(A) 
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observed  above,  it  now  follows  that  must  (also)  be  a  nonneg¬ 

ative  proper  combinatorial  basis  for  E(A )  satisfying  the  requirements  of  (a) 
with  respect  to  1Z(A).  Next  using  the  fact  that  Zat  -  ZT  >  0,  we  can 
follow  arguments  of  a  similar  spirit  to  those  used  in  showing  (a),  to  exhibit 
that  E(AT )  has  a  nonnegative  proper  row  combinatorial  basis  uM, . . . , 
satisfying  the  requirements  of  (b). 

We  come  now  to  the  proof  of  the  sufficiency  of  conditions  (i)  and  (ii). 
Let  uf1), . . .  and  uf1), . .  be  nonnegative  proper  column  and  row 

combinatorial  bases  for  E(A)  and  E{AT )  satisfying,  respectively,  the  re¬ 
quirements  of  (a)  and  (b).  First  note  that  for  any  i  ^  otj,cek,  where  i  6  (p) 
and  j,k  €  (m),  the  subvectors  and  cannot  be  concurrently 

nonzero  vectors.  For  if  that  were  possible,  then,  by  (ii),  i  is  must  be  a  non¬ 
singular  vertex  such  that,  in  ft(/l),  07..  y  i  y  otj,  which  is  not  compatible 
with  condition  (i).  This  shows  that  the  sets  of  vectors  {«(*), . .  .  and 

(tjf1), . . . ,  v(m)}  are  bi-orthogonal.  In  particular,  we  can  a  priori  scale  the 
vectors  in  these  sets  so  that  they  become  bi-orthonormal,  that  is, 

(vk)TxS ')  =  6k,h  k,l  6  (m).  (4.1) 

Define  the  matrix 

m 

y  =  Y^uUKvii))T-  (4.2) 

j=i 

Then,  as  can  be  readily  ascertained  using  (4.1),  Y  >  0,  Y  =  Y2,  and  the 
columnspace  of  Y  equals,  precisely,  E(A).  Moreover,  as  (v^)T  =  (yW)T Z 
for  all  j  6  (m),  we  see  that  Yw  =  0  for  any  generalized  eigenvector  of  A  cor¬ 
responding  to  an  eigenvalue  other  than  0.  Whence  Z  —  Y .  This  completes 
the  proof.  □ 

Our  theorem  shows  that  for  ZA  to  be  nonnegative  it  is  necessary  that 
E(A)  and  E(AT )  have,  respectively,  nonnegative  column  and  row  combi¬ 
natorial  bases  satisfying  its  requirements.  Therefore,  for  interest’s  sake  we 
provide  here  an  example  of  a  minus  M-matrix  A  such  that  E(A)  has  a  non- 
negative  proper  column  combinatorial  basis  satisfying  conditions  (ii)(a)  of 
the  theorem,  but  E(AT )  does  not  have  a  nonnegative  proper  combinato¬ 
rial  row  basis  satisfying  conditions  (ii)(b)  and  s^,  as  expected,  ZA  is  not 
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nonnegative.  Consider  therefore 


/-I  1  0  \ 

A\  =  I  1-10 

V  0  1  0/ 


Here  i/(Ai)  =  2  and  E(Ai)  =span{(  1  1  0)T,(0  0  1)T},  whereas 
E(Aj')  =span{(  1  1  0)  , (-1/2  0  l)r}.  A  computation  now  shows 


that 


ZAi  = 


1/2  1/2  0 
1/2  1/2  0 
-1/4  1/4  1  ) 


We  next  comment  that  Theorem  2  also  illustrates  that  the  conditions  for 
a  minus  M-matrix  A  to  possess  a  nonnegative  Perron  eigenprojection  are  not 
purely  combinatorial  in  the  sense  that  two  minus  M-matrices  can  have  the 
same  reduced  graph,  yet  one  possesses  a  nonnegative  Perron  eigenprojection 
while  the  other  one  does  not.  To  see  this  contrast  Aj  above  with  the  matrix 


M 


-i 

l 

1 


l  o  \ 

-1  0 

1  °/ 


Then  both  A\  and  A 2  ha,ve  the  same  reduced  graph,  but  whereas  Zax  is  not 
nonnegative,  we  find  that 


/  1/2  1/2  0  \ 

=  1/2  1/2  0  . 

V  0  0  1/ 

Additionally  we  mention  that  if  A  is  a  minus  M-matrix  for  which  Z  is  non¬ 
negative  and  in  71(A)  there  is  a  vertex  i  $  5(A)  with  access  to  a  vertex 
aj  G  5(A),  then  our  characterization  in  the  above  theorem  of  the  special 
basis  which  E(A )  must  possess  does  not,  a  fortiori,  tell  us  whether  the 
i-th  sub  vector  of  the  y-th  basis  vector,  namely,  is  zero  or  not.  That 
deoends  on  whether  Z:  a  is  zero  or  not.  Consider  the  followine;  examples: 
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A: ,  = 


0  0 

0  0 

0  0 

0  0 

-2  1 

1  -2 


/  0  0  0  \ 

A4  =  [  1/2  0  0  . 

\  1  1-1/ 

Observe  that  in  7£(/l3),  d(3, 1)  =  2,  and  that  in  7l(A4),  d( 3, 1)  =  2.  However, 
a  computation  shows  that 

1/2  1/2  0  0  0  0\ 

1/2  1/2  0  0  0  0 

0  01/21/200 
”  0  0  1/2  1/2  0  0 
0  0  1/2  1/2  0  0 

0  0  1/2  1/2  0  0/ 


— 


Thus  we  see  that  while  (^3)3,1  =  0,  (^4)3,1  >  0.  This  means  that  in  the 
basis  for  2?(.43)  provided  for  by  Theorem  2,  =  0,  whereas  in  the  basis 

for  E(A.{)  provided  by  the  theorem,  t/*)[3]  >  0. 

We  finally  remark  that  in  the  introduction  we  described  the  nonnegative 
basis  for  E{A)  which  exists  under  the  conditions  of  Theorem  2  as  “sparse”. 
This  is  because  when  contrasted  with  other  nonnegative  bases  which  exist 
for  E(A),  such  as  the  strongly  combinatorial  basis  given  in  (2.3),  we  see 
that,  at  least,  wherever  a  subvector  of  a  vector  in  the  latter  basis  is  positive 
because  of  an  access  from  a  singular  class  other  than  the  one  whose  index 
equals  to  the  index  of  the  basis  vector,  the  corresponding  subvector  of  the 
basis  vector  with  same  index  provided  under  the  conditions  of  the  theorem 
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5  ASYMPTOTIC  BEHAVIOR  OF  POWERS  OF 
A  NONNEGATIVE  MATRIX 


In  this  section  we  apply  our  Theorem  1  to  obtain  a  simple  proof  of  Theorem 
5.10  of  [5]  (see  also  Theorem  9.8  of  [12J.)  The  proof  given  in  [5]  depends 
on  analytical  results  concerning  convergent  and  summable  series.  We  shall 
work  with  an  n  x  n  nonnegative  matrix  P  which  we  shall  assume  to  be, 
without  loss  of  generality,  in  block  lower  triangular  Frobenius  normal  form, 
viz., 


(PiA  0  ...  0  \ 

I  P2,l  P2,2  0  I 


(5.1) 


\Pp,i 


Pp,p' 


where  the  diagonal  blocks  are  irreducible  or  the  lxl  null  matrix. 


We  shall  use  here  the  terminology  and  notations  of  [12,  Section  9]  with 
the  exceptions  as  noted  below.  We  call  a  vertex  i  of  72(P)  a  X-vertexU X  is  an 
eigenvalue  of  Pti;.  Now  let  i  and  j  be  vertices  of  7 2(P).  We  denote  the  spec¬ 
tral  radius  of  P{i,j)  by  s(i,  j)  and  we  put  A  =  P{i,j}  -  s(i,  j)/*  The  local 
distance  6(i,j )  is  defined  to  be  the  (singular)  distance  d(i,j)  in  A.  In  other 
words,  6(i,j)  is  the  maximal  number  of  s(i,  j)-vertices  in  any  path  from  i 
to  j  in  72(P).  (Note  that  6(i,j)  as  defined  in  this  paper  equals  d(i,j)—  1  as 
defined  in  [12]  or  [5].)  For  convenience  we  shall  let  Q  =  P{i,  j}.  For  an 
eigenvalue  A  of  Q  we  shall  denote  by  Z^(Q),  0  <  k  <  v\ (Q)  -  1,  the  k-th 
principal  component  of  Q  corresponding  to  A,  where  v\(Q)  is  the  index  of 
A  as  an  eigenvalue  of  Q ,  that  is,  its  multiplicity  in  the  minimal  polynomial. 
(Recall,  c.f.,  [9,  p.314]  that  Z^(Q)  =  AkZ^(Q),  0  <  k  <  v\(Q)  -  1,  where 
Z^{Q)  is  the  eigenprojection  of  Q  corresponding  to  A.) 

If  P  is  an  irreducible  nonnegative  matrix,  the  cycle  index  c(P)  of  P  is 
the  greatest  common  divisor  (g.c.d.)  of  the  lengths  of  all  simple  cycles  in  the 
directed  graph  of  P  (which  is  defined  in  [12]  and  many  other  papers.)  It  is 
well  known  that  c(P)  equals  the  number  of  eigenvalues  of  P  on  the  spectral 
circle  of  P.  Now  suppose  again  that  P  is  possibly  a  reducible  nonnegative 
matrix  in  Frobenius  normal  form  given  in  (5.1).  Let  i  and  j  be  vertices  in 
7 l(P).  We  now  define  the  quantity  g(i,j).  If  s(i,j)  =  0,  we  put  g{i,j)  =  1. 
Suppose  that  s(i,j)  >  0.  Let  II(i,j)  be  the  set  of  all  paths  from  i  to  j  in 
7 Z(P)  such  that  each  path  contains  6(i,j)  vertices  that  are  s(i,  j)-vertices. 
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Then  for  each  path  77  in  II(t,j),  let  <7(77)  be  the  g.c.d.  of  all  cycle  indices 
c(P*|fc),  where  k  ranges  over  all  s(i,  j)-vertices  of  77.  In  this  case  we  let  g(i,  j) 
be  the  least  common  multiple  of  all  <7(77)  for  paths  77  e  n(i,j).  Suppose 
s  =  s(i,  j)  >  0  and  put  g  =  g(i,j).  Define  now  the  smoothing  matrix 

=  p(P ), 


where  p(z )  is  the  polynomial 

p(z)  =  (l  +  2r/s  +  ...  +  (z/s)s_1)/(7.  (5.2) 

Thus  if  $  =  1,  the  smoothing  matrix  M(i,j )  as  defined  in  this  paper  equals 
M(iyj)/g  as  defined  in  [12].  (We  take  here  the  opportunity  to  point  out  two 
misstated  definitions  in  [12].  Definition  (9.7)  of  [12]  is  incorrect,  P  there 
should  be  replaced  by  P/s.  Also,  the  defintion  of  g(i,j)  given  in  [12]  has  an 
omission.) 

We  first  prove  the  following  lemma: 

LEMMA  4  Let  P  be  a  nonnegative  matrix  given  in  (5.1).  Let  i,j  £  (p) 
and  let  Q  =  P{i,j).  Suppose  that  6  =  S(i,j),  s  =  s(i,j),  and  g  =  g(i,j). 
Suppose  that  A  =  iu  6  a(Q) ,  where  |w|  =  1.  Then: 

(i)  (Q)  <6  =  u3(Q) 

and 


(ii)  if  v\(Q)  =  vs(Q),thenu9  =  1. 

Proof:  (i)  Let  C  be  a  path  from  i  to  j  in  7 2(P)  which  contains  a  maximal 
number  of  A-vertices.  Let  q  and  p  be,  respectively,  the  number  of  A-vertices 
and  s-vertices  on  £•  We  shall  show  that 

»\(Q)  <  q  <P  <  Vs{Q)  =  S. 

The  first  of  these  inequalities  is  immediate  from  Theorem  (5.9)  of  [7],  see 
also  [4].  To  prove  the  second  inequality,  we  note  that  it  follows  from  the 
Perron-Frobenius  theory  applied  to  the  irreducible  blocks  of  Q  that  every 
A-vertex  of  TZ(Q)  is  also  an  s-vertex  since  s  =  p(Q).  The  third  inequality 
is  immediate,  since  by  definition  6  is  the  maximal  number  of  s-vcrtices  on 
any  path  from  i  to  j  in  1Z(Q).  The  last  equality  follows  from  the  Rothblum 
index  theorem. 
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(ii)  By  the  Perron-Frobenius  theory,  w  is  a  root  of  unity,  say  it  is  the 
k- th  root  of  1.  Since  p  =  q,  every  5-vertex  on  £  is  also  a  A-vertex.  Hence 
(by  Perron-Frobenius)  the  cycle  index  of  every  vertex  on  £  is  a  multiple  of 
k.  It  follows  that  the  g.c.d.  of  the  cycle  indices  of  every  5-vertex  on  £  is  a 
multiple  of  k.  But,  since  p  =  q,  we  have  that  £  £  II(i,  j)  (which  was  defined 
previously  in  this  section).  Hence,  by  the  definition  of  g,  it  follows  that  g  is 
a  multiple  of  fc,  and  we  are  done.  □ 

We  are  now  ready  to  use  the  results  of  Section  3  to  re-prove  Theorem 
5.10  of  [5]. 


THEOREM  3  (Friedland  and  Schneider  [5,  Theorem  5.10])  Let  P  be  an 
nxn  nonnegative  matrix  in  given  in  (5.1).  Let  i,j  6  (p)  and  let  Q  =  P{i,j). 
Set  6  =  S(i,j),  s  =  s(i,j),  g  =  g(i,j),  and  let  M  =  M(i,j).  Suppose  that  i 
has  access  to  j  in  7£(/l).  Then  (Z(5-1)(Q)),-j  ^  0  an<l  if  s  >  0,  then 


(MPm)i,j  = 


+  o(m' 


5-1  em-5+l 


).  (5.3) 


Remark:  Note  that  the  matrices  Q ,  Z(Q),  and  N  used  in  the  statement 
and  proof  arc  block  matrices  with  fewer  blocks  in  their  bloc!;  dimension  than 
P.  We  shall  retain  the  block  row  and  column  indexing  of  P.  Thus,  for 
example,  Qjj  is  the  upper  left  hand  block  of  Q  and  Qjt;  is  the  lower  right 

hand  block  of  Q.  In  particular,  (zi5-1^(C2)),-(j  =  (ZsS~^(P))i,j- 


Proof:  By  Theorem  1,  (ZsS~l\Q))i,j  >  0.  Suppose  that  s  >  0.  We 
first  assume  that  s  =  1  =  p{P).  Let  N  =  p(Q),  where  p(z)  is  given  in  (5.2). 
By  Lemma  1,  (MPm),-j  =  (. NQ =  (Pm{Q))i,j,  where,  since  s  =  1, 

Pm(z)  =  zmp(z)  =  (zm  +  zm+1  +  ...+  zm+°-l)/g,rn  =  0,1,.... 


well  know  (cf.  Lancaster  and  Tismenetski  (9,  p.3M])  that 

*<r  =  E  E 

Aea(Q)  r=0 

Since 

=  M(?)A"'r  +  (mr  ^A”+,'r+---+(m+r  _1)A 


(5.4) 


m-r+g- 

(5.5) 
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it  follows  that 


"aMI-I  (r)/Av  /% 

*<r  =  E  E  &#2lr,«)+«d)- 


(5.6) 


A<=»(«),  |A|=1  ™o 


Since  by  Lemma  4(i),  we  have  that  v\(Q)  <  S  when  |Aj  =  1,  it  follows  from 
(5.5),  (5.6),  and  Lemma  4  that 


NQm  = 


(1  +  A  +  ...+  A 9-i)md-ixm-s+iz(s  D(g) 

1/  - ^TTv - +°(m  )• 

A6a(Q),  |A|=1  ^ 

(5-7) 

Suppose  that  A  €  a(Q),  |A|  =  1,  but  A  ^  1.  Then  either  X3  =  1,  in  which 
case  1  +  A  +  . . .  +  A3-1  =  0,  or  (Z^~1\Q))i,j  =  0  by  Lemma  4.  We  now 
obtain  that 


(NQm) 


.  _  l-'\Q))u  .  s.u 


(5-1)! 


(5.8) 


To  obtain  the  general  case  of  s  >  0  from  (5.8),  we  note  that  pm(sQ )  = 
smpm(Q)  and  that  Zj-l(sQ)  =  ss~lzf~l\Q).  □ 


We  observe  that  if  s(i,j)  =  0  for  a  nonnegative  matrix  P  given  by 
(5.1),  then  P{i,j }  is  strictly  lower  triangular  and  hence  nilpotent.  Thus 
(Pmki  =  o  for  sufficiently  large  m. 
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